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A generating coordinate is introduced into the exchange-correlation functional of density- 
functional theory (DFT). The many-body wave function is represented as a superposition of Kohn- 
Sham (KS) Slater determinants arising from different values of the generating coordinate. This 
superposition is used to variationally calculate many-body energies and wave functions from solu- 
tions of the KS equation of DFT. The method works for ground and excited states, and does not 
depend on identifying the KS orbitals and energies with physical ones. Numerical application to the 
Ifelium isoelectronic series illustrates the method's viability and potential. 
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Density-functional theory (DFT) is routinely used in 
quantum chemistry and physics for the calculation of, 
e.g., ground-state energies and charge distributions. In 
a typical application of DFT one solves the Kohn-Sham 
(KS) single-particle equation, with a suitable choice for 
the exchange-correlation functional, and uses the result- 
ing single-particle orbitals to construct the ground-state 
particle density. From this density other observables, 
such as the ground-state energy, can then be calculated. 
It is well known that the KS Slater determinant is not 
in itself an approximation to the many-body wave func- 
tion but only a device for reproducing the correct density. 
Similarly, the single-particle eigenvalues of the KS equa- 
tion do not represent the energy spectrum of the system 
under study, and their differences are not excitation en- 
ergies. Although excited-state energies and many-body 
wave functions for both ground and excited states are, in 
principle, functionals of the ground-state density,^ they 
are hard to extract from DFT. For excited-state energies 
a number of viable methods has been suggested, ^'^'"^'^ but 
these are more complicated than conventional ground- 
state DFT, and not yet as widely used. Many-body wave 
functions cannot be calculated at all from standard KS- 
type calculations. 

In the present paper a novel approach based on DFT is 
proposed, which allows to (i) systematically improve the 
accuracy of ground-state energies obtained in a standard 
KS calculation, (ii) calculate excited-state energies, and 
(iii) obtain variational approximations for the many-body 
wave functions corresponding to these energies. The de- 
vice which makes this possible is the GrifRn-Hill- Wheeler 
(GHW) variational method,®'^ or generator-coordinate 
method (GCM), which is now briefly described. In the 
GHW approach one writes the trial wave function 4" as 
an integral transformation of a generating wave function 
4>, according to 



(1) 



where a is a generating coordinate and the variational 
principle is applied to the kernel of the integral, f{a). 
The nature of 5* and $ (i.e., whether they are many- 
body or single-body functions) does not matter at this 
stage, and their coordinate arguments have thus been 
suppressed. Calculation of the expectation value of the 
system's Hamiltonian H with followed by variation of 
the result with respect to f{a), yields the so-called GHW 
equation^ 



da' [K{a, a') - ES{a, a')] /(a') = 0. 



(2) 



Here K{a,a') = ($(a)|i/ |$(a')) and S{a,a') = 
(<i>(a)|<i>(a')) are matrix elements of the Hamiltonian 
with the generator functions $(«), and the overlap of 
these functions, respectively. Solution of the GHW equa- 
tion (2) yields the energies E, which are variational ap- 
proximations for the eigenvalues of H, and the function 
f{a) which through (1) determines the corresponding 
wave functions. The full set of solutions of the eigenvalue 
problem (2) yields thus the spectrum and wave functions 
of Hamiltonian H. 

Originally this method arose in nuclear physics, ^'^ 
where 4' was taken to be the nuclear many-body wave 
function, and the generating coordinate a was inter- 
preted as a 'deformation parameter' describing collec- 
tive oscillations of the nucleus. The generator functions 
<i>(a) were obtained from solving an auxiliary, simplified, 
Schrodinger equation in which the nuclear potential v(r) 
was replaced by a deformed potential Va (r) , with a char- 
acterizing the degree of deformation. Outside nuclear 
physics the GHW approach has been applied to vari- 
ous model Hamiltonians,® the electron gas,^ and molec- 
ular electronic-structure calculations, ^"^'^^'^^ exploring a 
variety of different choices for the generator functions 
<& and the generating coordinate a. More recently, it 
has been used to generate optimized basis functions for 
Hartree-Fock calculations, starting from a set of simple 
trial functions. "'^^^^'^'"'^^ In this type of application a is iden- 
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tified with the basis function exponent C, and ^' and $ 
are single-particle functions. 

In the present contribution we go back to many-body 
wave functions. The basic idea is to identify the auxil- 
iary Schrodinger equation with the KS equation of DFT, 
the deformation potential producing the family of genera- 
tor functions with the KS potential, and the deformation 
parameter (or generating coordinate) with a parameter 
in the exchange-correlation (xc) functional. Many func- 
tionals, such as the Xa and B3-LYP approximations, 
naturally contain such a parameter, but it can always 
be introduced by hand in any functional. The present 
proposal is thus to write the many-body wave function 
as 

*(ri,...rjv) = y"da/(a)$^^(a;ri,...rjv), (3) 

where $^'^(a) is the Slater determinant obtained from 
a KS calculation with xc potential Vxc,a- Note that the 
deformation parameter a is neither an adjustable param- 
eter fixed by comparison with experiment nor a varia- 
tional parameter for a given form of the trial wave func- 
tion, but rather a generating coordinate that accounts 
for collective behaviour in not described by the single- 
particle coordinates ri . . .r^r in the Slater determinant 
[a; ri . . . yn)- Its physical origin is in the deforma- 
tions of the single-body potential Vxc,a, which simulate 
the collective degrees of freedom of the interacting many- 
body system.^ 

The idea proposed here is thus to use DFT poten- 
tials and orbitals as input for a GHW calculation, and 
the GHW Eq. (2) with (3) to gain direct variational 
access to many-body energies and wave functions. As 
a first viability test of this scheme it is now applied 
to the Helium isoelectronic series. For a general two- 
electron atom in a closed-shell configuration the ker- 
nels K and S are evaluated easily. $^'^(a) is a 2 x 2 
Slater determinant formed with one doubly occupied or- 
bital. After performing the sum over spins one is left 
with only spatial integrals, which need to be calculated 
numerically, and the kernels take the form K{a^ a') = 
2{a\a'){a\t + ^a') + {aa\u\a' a') , Sld.A S{a,a') ~ (aja')^, 
where \a) stands for the KS orbital ipa{r), and t and ij 
are the single-particle kinetic and potential energy oper- 
ators, u = l/|r — r'|. Note that the diagonal element 
[a = a') of the kernel K is simply the energy expres- 
sion one obtains in a restricted closed-shell Hartree-Fock 
calculation^^ for a two-electron atom. A similar formal 
connection to Hartree-Fock theory will always hold, since 
K is the matrix element of the Hamiltonian between sin- 
gle Slater determinants. The single-particle orbital Lp^ (r) 
is now obtained from a self-consistent KS calculation. 
For simplicity (and to illustrate an important point be- 
low) I choose as xc functional the Xa approximation and 
take as generating coordinate the parameter a present in 
that functional. (The fact that the GHW generating co- 
ordinate and the Xa coefficient arc both traditionally 
called a is a coincidence.) The GHW integral equation 



TABLE I: Ground-state energies of the Helium isoelectronic 
series. First row: energies obtained from solving the GHW 
equation, using the Xa functional to obtain the generator 
orbitals, as described in the main text. Second row: exact 
(nonrelativistic) energies from Ref. 18. Third row: percentual 
deviation of the GHW energy from the exact one. All energies 
are in atomic units. For comparison: the LDA result for He 
is = —2. 835a. w. and deviates from the exact value by 

2.4%. 







rriexact 


% deviation 


He 


2.871 


2.904 


1.1 


Li+ 


7.244 


7.280 


0.49 


Be^+ 


13.62 


13.66 


0.29 




22.00 


22.03 


0.14 




32.37 


32.41 


0.12 




44.75 


44.78 


0.067 


0"+ 


59.12 


59.16 


0.068 




75.50 


75.53 


0.040 



(2) is then solved by discretization. A suitable mesh for 
discretizing the integral over a is {0, 0.5,1, 1.5, 2}. Five 
mesh points may appear surprisingly few, but empiri- 
cally it was found that larger or denser meshes did not 
significantly change the results. The calculation thus 
consists in the following three steps: (i) performing KS 
Xa calculations at the prescribed values of a; (ii) using 
the resulting orbitals to evaluate the kernels K and 5*, 
according to the above equations; and (iii) discretizing 
the integral equation (2) on this set of a's and solving 
it by standard matrix algebra. In the remainder of this 
paper some representative results obtained in this way 
are presented. 

Table I compares the lowest eigenvalue of Eq. (2) (i.e., 
our approximation to the many-body ground-state en- 
ergy) with reference data obtained from numerically ex- 
act wave functions. Considering the simple discretiza- 
tion scheme and generating functional {Xa), the quan- 
titative agreement achieved is rather surprising. The 
largest deviation from the reference data is found for He- 
lium, and is only about 1%. Interestingly, this agreement 
has been obtained by starting out with an xc functional 
that on its own yields significantly worse energies: the 
ground-state energy of He calculated from the usual KS 
scheme, employing the Xa. functional with the above val- 
ues of a, is found to be E'f^=^ = -1.952a.M., E^=^-^ = 
-2.515a.u., E^='^ = -3.170a.M., E^='^-^ = -3.915a.u., 
and Eq^"^ = — 4.749a.u., respectively, which are all sig- 
nificantly off the true He ground-state energy. Out of 
these unphysical energies the GHW optimization gener- 
ates a ground-state energy of Eq^^~'^°' = — 2.871o.u., 
which deviates only by about one percent from the exact 
reference value (see Table I). Since this is a variational 
calculation it is, of course, rather natural that the GHW 
energy is better than that obtained with each generating 
functional (i.e., each value of a) individually. GHW opti- 
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mization can thus be used to systematically improve on 
results obtained from a given input density functional, 
which need not be very good on its own. On the other 
hand, the Xa GHW value — 2.871a.u., found above, is 
also closer to the exact result^^ j^exact _ _2.904a.t(. 
than the ones obtained with Hartree-Fock^* {E^^ = 
-2.8620.U.) or LDA (i^^^^ ^ -2.835a.u.)- This im- 
provement may appear surprising because the Xa func- 
tional in itself can be interpreted as a rather simple ap- 
proximation to both Hartree-Fock and LDA, but is sim- 
ply explained by noting that in the present context the 
Xa functional (or any other generator functional from 
DFT that could be used instead) only serves as a conve- 
nient way to build a family of continuously parametrized 
generator determinants $(a), and is not used directly to 
obtain the desired results. 

In principle, the N eigenvalues of the GHW equa- 
tion found by discretizing it as 'aw N x N matrix equa- 
tion provide the N lowest-lying energies of the original 
Hamiltonian, but the quality of the resulting energies 
will depend on the nature of the generator functions. 
The present calculation takes ground-state KS determi- 
nants to form the generator function and thus does not 
directly aim at excited states. However, from the second- 
lowest eigenvalue of the GHW equation one obtains an 
estimate for the energy of the lowest excited state with 
same symmetry as the ground state. For He this is 
the 2^5*0 para state, with an energy of — 2.146o.u.. The 
value found in the above ground-state GHW calculation 
is E^"^-^°' = -1.7880.U. and deviates from this by 
16.7%. A ground-state GHW-DFT calculation thus al- 
lows one to obtain estimates of the energies of excited 
states with the same symmetry as the ground state. It 
is thus complementary to methods based on minimiza- 
tion in symmetry subspaces, which give access to ex- 
cited states of symmetry different from the ground state. 
Calculations optimized for specific target excited states, 
using excited-state KS Slater determinants as generator 
functions instead of ground-state ones, are expected to 
provide better results for the corresponding excitation 
energies. 

Another interesting application of GHW variational 
optimization in DFT is the calculation of many-body 
wave functions. Within the GHW scheme the many-body 
wave function ^ is known in terms of the family of gen- 
erator functions $(a) and the weight function /(a), by 
means of Eq. (1). After discretization, the weight func- 
tion reduces to a set of coefficients /(a), which can be 
read off directly from the components of the eigenvector 
corresponding to a given eigenvalue E of Eq. (2). These 
coefficients are thus automatically obtained together with 
the eigenvalues. As an explicit example, the (unnormal- 
ized) many-body wave function obtained together with 
the above result of E§"^ = -2.871a.w. for the He 
ground state is 

= -0.0529$^io + 0.276$^io.5 - 0.451$^ii 

-h0.770$^ii.5 - 0.354$^i2> (4) 



TABLE II: Expectation value of the operator r" for various 
values of n calculated with the Helium GHW many-body wave 
function (4), compared with results from a standard DFT 
calculation using the B88-LYP functional^^'^° (in a.u., and 
after normalization). 

n {r'^)GHW-Xa-DFT {t^) BSS-LY P-DFT 

-2 5.74 5.98 

-1 1.66 1.69 

1.00 1.00 

1 0.928 0.964 

2 1.16 1.26 



while the one obtained together with Eq^^ = 
—59. 12a. M. for the ground state is 

= -0.116$^fo + 0-539$fio.5 - 0.786$^fi 

-F0.279$^ii.5 + 0.0263$^i2, (5) 

where the (^^^{a) are the KS Slater determinants ob- 
tained in the generator Xa calculations for He and , 
respectively. Wave functions for excited states are ob- 
tained in just the same way from the components of the 
higher eigenvectors. 

In a variational calculation, such as the one performed 
here, energies are typically obtained with higher accu- 
racy than wave functions. However, the simple five-term 
expansions, given above, already suffice to obtain expec- 
tation values of observables that are comparable to those 
obtained with other methods. As an example, consider 
the expectation value (5'o|r"|^'o), for n ranging from —2 
to +2. Results for the He atom, calculated with the 
GHW wave function (4) are listed in Table II. In spite 
of the simplicity of the Xa functional and the inferior- 
ity of variational wave functions to variational energies, 
the expectation values obtained with the five-term GHW- 
optimized Xa wave function are found to be close to 
those obtained from a much more sophisticated density 
functional. 

Expressions (4) and (5) for '^q also allow us to see 
clearly the difference between the present scheme and 
configuration interaction (CI). In CI ^'o is written as a 
linear combination of Slater determinants that all stem 
from the same HF calculation. Individual determinants 
differ by systematically substituting occupied single- 
particle orbitals by unoccupied ones. In the present 
scheme is also a linear combination of Slater deter- 
minants, but each determinant comes from a different 
KS calculation, the amount and nature of the difference 
being specified by the deformation coordinate a. The 
determinants in the present expansion thus stem from 
Hamiltonians with different potentials. This implies that 
each of them can be interpreted as an effective rcsumma- 
tion of a large number of Cl-type determinants arising 
from a fixed Hamiltonian. Another consequence is that 
the individual determinants in the GHW expansion are 
not mutually orthogonal. Interestingly, representation of 
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many-body wave functions in terms of non-orthogonal 
determinants has previously been shown, in the context 
of the resonating Hartrec-Fock method, to be an efficient 
way to account for strong Coulomb correlations not read- 
ily accounted for by traditional expansions. ^^'^^'^^ 

In summary, the present paper proposes to combine 
two many-body methods (DFT and GHW), each of which 
is successful in its field of origin (electronic-structure the- 
ory and nuclear physics, respectively), but which had not 
previously been brought to work together. At the heart 
of the present proposal is Eq. (3), which expresses the 
many-body wave function as a weighted superposition 
of Kohn-Sham determinants, each arising from a differ- 
ently deformed exchange-correlation potential. As a first 
viability test the method has been applied to the He 
isoelectronic series. A full judgement of the powers of 



GHW-DFT must await systematic tests for a wide va- 
riety of physical systems and generator functionals, but 
the present initial exploration shows that (i) ground-state 
energies can be obtained that are considerably better 
than those calculated from both the generator functional 
and more sophisticated methods approximated by it, (ii) 
many-body wave functions are obtained with almost no 
additional numerical effort, and yield expectation values 
that are close to those given by other methods, and (iii) 
access to excited state energies and wave functions is, in 
principle, possible, too. 

Acknowledgments This work was supported by 
FAPESP. I thank L. N. Oliveira, M. Teter and A. Niklas- 
son for useful remarks on an earlier version of this paper, 
and M. Trsic and A. B. F. da Silva for hospitality at the 
IQSC. 



* Electronic address: capelle@iqsc.sc.usp.br 

^ P. Hohonborg and W. Kohn, Phys. Rev. 136, B864 (1964). 
^ E. K. U. Gross, L. N. Oliveira, and W. Kohn, Phys. Rev. 

A 37, 2805; ihid 2809; ihid 2821 (1988). 
^ M. Petersilka, U. J. Gossmann and E. K. U. Gross, Phys. 

Rev. Lett. 76, 1212 (1996). 
" A. Nagy, Chom. Phys. Lett. 296, 489 (1998). 
^ A. Gorling, Phys. Rev A 59, 3359 (1999). 

J. J. Griffin and J. A. Wheeler, Phys. Rev. 108, 311 (1957). 

D. L. Hih and J. A. Wheeler, Phys. Rev. 89, 1106 (1953). 

* P. Chattopadhyay, R. M. Dreizler, M. Trsic, and M. Fink, 
Z. Phys. A 285, 7 (1978). 

® B. Johansson and J. da Providencia, Physica 94 B, 152 
(1978). 

B. Laskowski, J. Diamond, A. Walch, and B. Hudson, J. 
Chcni. Phys. 69, 5222 (1976). 

B. Laskowski in Quantum Science: Methods and Structure, 
J.-L. Calais, O. Goscinski, J. Lindberg, and Y. Ohrn, eds. 

(Plenum, Now York, 1976). 

E. Dcumons, L. Lathouwcrs. P. van Louvon, and Y. Ohrn, 
Int. J. Quantum Chcm. Suppl. 18, 339 (1984). 

^■^ J. R. Mohallem, R. M. Dreizler, and M. Trsic, Int. J. Quan- 
tum Chem. Symp. 20, 45 (1986). 

1* A. B. F. da Silva, H. F. M. da Costa, and M. Trsic, Mol. 
Phys. 68, 433 (1989). 

15 A. B. F. da Silva, G. L. MalU, and Y. Ishikawa, Chem. 



Phys. Lett. 203, 201 (1993). 

A. Szabo and N. S. Ostlund, Modem Quantum Chemistry 
(McGraw-Hill, New York, 1989). 

It may also be recalled in this context that in both single- 
particle applications of the GHW scheme"'^*'^^'^*'^'^ and 
model calculations* good results have been obtained with 
a comparably small number of discretization points. A ra- 
tionale for the success of simple discretization schemes has 
been proposed in Refs. 24,25. 

E. R. Davidson, S. A. Hagstrom, S. J. Chakravorty, 
V. M. Umar, and C. Froese Fischer, Phys. Rev. A 44, 
7071 (1991). 

A. D. Becke, Phys. Rev. A 38, 3098 (1988). 

C. Lee, W. Yang, and R. G. Parr, Phys. Rev. B 37, 785 

(1988). 

H. Fukutome, Prog. Theor. Phys. 80, 417 (1988), ibid 81, 
342 (1989). 

N. Tomita, S. Tenno, and Y. Tanimura, Chem. Phys. Lett. 
263, 687 (1996). 

S. Tenno, Theor. Chem. Acc. 98, 182 (1997). 

J. Brocckhove and E. Deumens, Z. Phys. A 292, 243 

(1979). 

F. Arickx, J. Broeckhove, E. Deumens, and P. van Leuven, 
J. Comp. Phys. 39, 272 (1981). 



